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SUBSONIC FLOW OVER THIN OBLIQUE AIRFOILS AT ZERO LIFT 

By Robert T. Jones 


SUMMARY 

A previous report gave calculations for the pressure distri- 
bution over thin oblique airfoils at supersonic speed. The 
present report extends the calculations to subsonic speeds. 

It is found that the flows again can be obtained by the super- 
position of elementary conical flow fields. In the case of the 
swe pi-back wing the pressure distributions remain Qualitatively 
similar at subsonic and supersonic speeds. Thus a distribution 
similar to the Ackeret type of distribution appears on the root 
sections of the swept-back wing at M=0. The resulting posi- 
tive pressure drag on the root section is balanced by negative 
drags on outboard sections. 

INTRODUCTION 

So far as is known, attempts to extend airfoil pressure- 
distribution calculations to three-dimensional flow hare 
been confined to cases of thin lifting surfaces. It has gen- 
erally been assumed that the component of the pressure 
distribution arising from the thickness of the airfoil will be 
but little affected by the finite span, or aspect ratio, of the 
wing. This supposition is borne out by the known incom- 
pressible-flow solutions for flat ellipsoids. These solutions 
show that the usual variations of aspect ratio produce small 
effects. 

Compressible-flow theory shows, however, that the effects 
of plan form become more pronounced at higher speeds. 
The theory indicates a progressive reduction of the equiva- 
lent aspect ratio as the Mach number approaches 1.0. 
Hence at these speeds the three-dimensional character of the 
flow can no longer be neglected. Of particular interest are. 
the deviations from two-dimensional flow near the root 
sections of a swept-back wing, since the adverse effects of 
compressibility may arise first in this region. 

In the present report three-dimensional flows are obtained 
from a distribution of “pressure sources” in the chord plane 
of the airfoil. The shapes thus obtained are symmetrical 
airfoils at zero lift. The calculations are simplified by 
considering airfoils composed of conical or cylindrical surfaces. 
In these cases the sources can be arranged into lines of uniform 
strength following the generators of the surface. The relation 
between the strengths of the line sources and the shape of 
the airfoil is the same as in reference 1; that is, each line 
source produces a deflection of the streamlines crossing over 
the source. The pressure field of the line source again can be 
represented by systems of straight rays of equal pressure 
(isobars) radiating from the ends of. the line source. 


In general, the present development follows closely that of 
reference 1 and the reader should consult that report for 
additional details of the method. The solutions are given 
explicitly for Af=0 but are extended to other Mach numbers 
by the well-known Prandtl transformation. 

THE OBLIQUE LINE SOURCE 

It is well known that an individual velocity component of a 
potential flow will satisfy the same differential equation as the 
potential. In the approximation of the thin-airfoil theory 
the pressure depends only on the individual component u, 
that is, 

2 V 


while the slope of the surface depends only on the individual 
component w, that is 


dz_w 

Tx~V 


( 2 ) 


(See appendix for symbols.) Hence in the thin-airfoil theory 
it is often more convenient to deal directly with the velocities 
u and w as solutions of Laplace’s equation than to derive 
these components from a velocity potential <p. 

Since u is proportional to the .pressure, a solution of 
Laplace’s equation can represent directly the pressure distri- 
bution, hence the term “pressure potential.” In this 
terminology, the fundamental solution 

u =f = fA+y r +? (3 > 

represents a point source of pressure rather than a point 
source of fluid. 

To get the effect of a row of sources, or a line source, along 
the x axis between the points a and b, it is necessary to inte- 
grate equation (3): 

_ P 

U V (*— !)*+#*+ z 2 

=sinh- ! ;Jq=i- sinh- 1 ' ...(4) 

The pressure field of the finite line source thus consists of the 
sum of two conical pressure fields radiating from the ends of 
the line source. (See fig. 1 .) In the supersonic case (refer- 
ence I), the radial isobars forming the conical field were 
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confined to the downstream Mach cone. Hetre, however, 
the isobars extend over the whole space. 1 

If the direction of flight is along the axis of the source 
(x axis), the flow will satisfy the boundary condition for a 
body of revolution. However, if the line source is tinned 
out to a position oblique to the stream, the boundary shape 
will be distorted and, if the angle of obliquity is large enough 
to place the line source well outside the diameter .of the 
original body, the figure formed will be an oblique wedge. 
The nose angle of the wedge is formed .where the streamlines 
of the main flow cross the line source. - 
At supersonic speeds the expression for the oblique line 
source was obtained by applying an equivalent of the Lorentz 
transformation, for which the wave equation is invariant. 
The equivalent transformation for Laplace’s equation is a 
rotation of the axes, given by 

x'=x-\-my 
y'=y—mx 
s! =z-yfl+m? 

1 The conical pressure field for either the suhsonlo or the supersonic lino source may be ob- 
tained directly from the general solutions of Laplace’s equations of sero degree in x, v, t given 
by W. F. Donkin. (See reference 2, page 867.) The general solution Is 

wfarew) 

The solution corresponding to the subsonic line source is 

_ _ ffW* 4 x 

Jog ■ . — / • 

s-f-yk’+fH-s 2 - Vy»+2» 

while the field for the supersonic source Is given by 

t+iz x 

ti“ — R, P, jog . — sas 1 ■ - ^ " 

r+^/r3— 1 


where m is the slope of the new axes relative to tho old. 
(Note that a change of scale is admitted for convenience.) 
The geometry of the pressure field relative to the line source 
is not altered in any way by this rotation and the isobars 
behave as though they were rigidly attached to the ends of 
the source. For a line source with one end at the origin, 
we have 


w=sinh _I 




"(5) 


This field is illustrated in figure 2 for the plane z— 0. As 
m— » oo the x and y axes interchange and there is obtained 


■ - (0) 

for a line source along y. 

The vertical velocity w near 2 = 0 , winch determines the 
shape of the boundary; may be found by integrating u with 
respect to x and then differentiating the resulting velocity 
potential with respect to z. 


d<f> 5 C z , 
u>=k— =5— u ax 
2)2 52 J-o. 


(7) 


Evaluation of this integral for the overlapping fields from 
two ends of a line source gives 


w=±2ic 


Vl-t-m 8 

m 


over the area of the xy plane behind the line source. 


( 8 ) 
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The figure formed by the streamlines crossing a line source 
is thus a wedge-shaped body having an oblique leading edge 
and extending indefinitely downstream. It is evident from 
equation (4) that the infinitely wide wedge cannot be treated 
in subsonic flow, since it creates an infinite pressure disturb- 
ance at all points. 

The slope of the wedge surface away from the chord plane 
is given by 


dz_w 

dx~V 


( 9 ) 


With this relation and equation (8) the pressure coefficient 
near the plane z — 0 may be expressed in terms of the slope 


Ap__ 1 to dz 

s * vr+m 5 



,x'-V 

~wr 


sinh -i 


I y'\ ) 


( 10 ) 


where \y'\ indicates the absolute magnitude of y' . Following 
the thin-airfoil theory, the pressure over the chord plane 
( 2 -» 0) is taken as the pressure over the actual airfoil surface. 


AIRFOILS BOUNDED BY PLANE SURFACES 

It was seen that the effect of a line source in the pressure 
field is to causo a deflection of the streamlines crossing the 
source. The deflection thus produced is. equal and opposite 
at points above and below the chord plane, so that the 
source spreads the streamlines apart. If the source is 
followed by a sink of equal strength, an equal opposite deflec- 
tion of the streamlines will occur as they cross over the sink. 
The figure formed by the streamlines near the plane z= 0 
will thus be a plate of uniform thickness with a beveled 
leading edge. 2 (See fig. 3.) 

• According to the thln-alrfoil theory the thickness of the flguie ends abruptly at the ends 
o£ the source lines. A more met consideration would be expected to show some rounding 
at the tips of the wedge as indicated In figure S. 
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Tlie pressure distribution over such a beveled edge may 
be obtained very simply by superimposing the pressures laid 
off on radial isobars originating from the four comers of the 
bevel. Figure 3 illustrates this process for a bevel having a 
square plan form. Only isobars from one tip are shown 
because of the symmetry of the figure. 

In figure 3, the line source and the line sink are parallel 
to the y axis, hence 


u=sinh - 


i y+1 

T*+i| 




sinlr 1 


y+1 

k~ If 


-j-sinh -1 


y — i 

\z-l[ 


( 11 ) 


It can be seen that if the aspect ratio of the figure is increased 
to a large value the ends of the line sources will be separated 
by a great distance and the isobars in intermediate regions 
will approach parallel straight lines, hence the flow field 
approaches a cylindrical or two-dimensional form. At the 
same time the arguments y±l/[a;±l| in equation (11) be- 
come i/±ri/x±l and ij takes on very large values so that 


■“■‘lffir * 1 '’ 82 


y±v 

(r+T[ 


and equation (11) is found to approach the Legendre func- 
tion Q 0l that is 


•w.=2 log 


x-1 


z+1 



4 Q 0 (x) 



(See reference 3, p. 110.) 

This expression when combined with equation (8) agrees 
with the two-dimensional potential function for the wedge, 
that is, 


— (u— iw) =ri Q 0 (x) ±2vi P 0 (x) (13) 


(See fig. 4.) 

The isobars at right angles to the axis of the line source are 
lines of zero pressure, hence the rays originating at the tip 
of a rectangular wing contribute nothing to the pressure 
distribution at this tip. The whole pressure distribution 
at one tip is thus obtained by considering only those isobars 
radiating from the opposite tip. It is evident that in the 
case of a long narrow rectangular wing the pressures at 
either tip will be approximately one-half the pressures over 
the middle portion of the wing. 

In case the wing is oblique the tip sections will no longer 
be at right angles to the axes of the source lines and the rays 
originating from the adjacent ends of the source lines will 
contribute to the pressure over the tip. It can be shown 
that this component of the tip pressure distribution is 
similar in form to the Ackeret type of distribution, that is, 
the pressure at any point of the surface is proportional to 
the slope of the surface at that point. 

Consider first the sloping surface formed by a pair of 
oblique source-sink lines. The tip section lies along the 
lines of constant pressure of magnitude proportional to 
sinh -1 lfm. Between the source and sink the pressures are 
additive, so that 


Ap 2 m az ... 1 
— = - r ~ — = -j- sinh 1 - - 
2 v VT+nr ax m 


(14) 


Ahead of or behind this section the pressures cancel. 

In ease of a curved airfoil surface the chord can be divided 
into elements composed of source-sink pairs, the strengths of 
which are proportional to the slope of the surface at the point 
in question. Each pair then contributes a pressure propor- 
tional to the local slope and contributes no pressure at other 
points. Hence, equation (14) applies when dzjdx is variable 
along the chord, - 




Figure 4.— Two-dimensional velocity functions for wedge and biconvex sections. 


The foregoing arguments of course apply only at- the tip 
section of the oblique wing. At some distance from the tip 
section the overlapping isobars radiating from the tip again 
produce a quasi-cylindrical pressure field as in the case of 
the rectangular wing. Thus the resultant pressure distribu- 
tion at either tip of a long oblique wing consists of two com- 
ponents, one given by equation (14) and of the Ackeret type 
while the other component is equal to one-half the normal 
two-dimensional pressure distribution associated with the 
airfoil section. 

Figure 5 shows the pressures over a beveled-edge profile 
having 45° sweepback. The pressure distribution over the 
root section is given by 


Ap 

—4 m dz\ 

2 . 

w Vl+ m * da: L 


£&>(:») -sinh' 1 ^ Po(z)J (15) 
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at a great distance from either root or tip by 
A v —4 m dz „ , . 

<L~ * VTW Tx Qo{x) 

and at the tip by 


Ap —2 

m dz f 

2 " r 

Vi+m*dz L 


Pian 



(16) 


( 17 ) 



Figcre 5.— Pressure distrlbattou over beveled edge vrith 4S° svreepback. 


To take account of the effect of compressibility we make 
use of the Prandtl transformation, increasing both the x 


dimensions and the pressure coefficients by the factor 


1 

-JT=W' 


Replacing m by -Jl—AP cot A, where A is the angle of sweep- 
back, equation (16) reduces to 


Ap —4 1 dz __ 

p/2(F cos A) 2 ir Yl — PTcoTA)* 5^coTa) 


Qo(x) (18) 


Thus, at a great distance from either root or tip, ihe pressures 
follow a variation indicated by the normal component of 
velocity V cos A. 

At the root section, a component representing the Ackeret 
type of pressure distribution is added to equation (18) . This 
component is 

4 1 dz 

ir 7 1— (M cos A)* d{x cos A) 

(I9) 


The factor sinhr 1 


shows a logarithmic in- 


Vl —M* cot A 

ffnity at M=1.0. Hence the pressure on the root section 
increases more rapidly with Mach number than do the pres- 
sures at other sections of the swept-back wing. Further- 
more, the shape of the pressure distribution over the root 
section approaches the Ackeret shape more closely as the 
Mach number approaches 1.0. As shown in reference 1, the 
pressure distribution on the root section is exactly this shape 
at supersonic speeds, that is, 


A p _4 1 dz 

2 cos* A x Vl— (itf cos A) 2 <2 ( x cos A) 

(20) 

Since sinh -1 -»cosh~ l for large values of the argument, the 
swept-back airfoil shows no discontinuity in the type of 
pressure distribution on passing through the speed of sound. 
It will be evident that similar reasoning can be applied to 
the tip sections. 

AIRFOIL OF BICONYEX SECTION 


The use of a finite number of sources and sinks results in 
airfoil sections composed of straight segments. Such sec- 
tions are undesirable, since they show infinite pressure peaks 
at the bends in the surface. Surfaces having continuous 
curvature require continuous distribution of sources and sinks 
alined with the generators of the surface. The simplest of 
these is the biconvex profile in which the upper and lower 
surfaces are parabolic arcs and have Constant curvature. 
Such a profile requires line sources of finite strength to form 
the desired angles of intersection of the arcs at the leading 
and trailing edges together with a uniform distribution of 
sinks along the chord plane between the two sources. 

The pressure field for a uniform sheet of line sources is 
obtained by integrating the field of a single line source in 
the x direction. This integral is 

J 5 ^- 1 #!*- sUr ‘ fr 

( 21 ) 

The integration for a source sheet is actually somewhat 
simpler if the interference of a bilaterally symmetrical ar- 
rangement of sources is considered simultaneously. The in- 
fluence of the symmetrical, or conjugate, arrangement is 
obtained by substituting — m for m in equation (21). De- 
noting x—my by x' and y+mx by y' we have 

^ («+«) = J (sinh- 1 sinh -1 dx 

=b (y' sinii_I m ~ y ' 8inh_1 ifr) 22) 

To obtain a complete swept-back wing it is necessary to 
add a number of component pressure fields as explained in 
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reference 1. For an infinite swept-back wing with leading 
and trailing edges at y'=-\-m and —m, respectively, on one 
side, and at ~y'= +m and — m, respectively, on the other side, 
there is obtained 


mm w / Binh -1 

2 IT Vl+m* \c/m«L^\ 


x'+l 


(23) 


\y'—m\ 

9 tah "‘lFT^) + l( 8 blh ''ir=ST _ ■ 

where ( - ) is the thickness-chord ratio of the biconvex 

profile. The terms Q x represent the pressure distribu- 
tion on the biconvex airfoil in two-dimensional flow. The 
appearance of these terms is the result of the assumption 
that the tips are removed to a great distance. 

At the root section (y=Q) equation (23) reduces to 


Ap 
2 ‘ 


■MsGL [ 4ft< *>- 4 8inl, “ s-™] (24> 


Figure 6 shows pressure distributions at various stations 
along tiie span for a biconvex wing with 60° sweepback. 
The curves assume the two-dimensional form at a relatively 

short distance (y^~c) from the root section, and similar 
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Figure 6 —Pressure distribution at various spar ’vise stations on swept-back wing, A -60°, 

At— 0. 


behavior is to be expected near the tips. Hence the assump- 
tion of. infinite aspect ratio should apply very nearly at any 
section situated more than one-half chord length from either 
root or tip. 

Figure 7 shows the effect of Mach number on the pressures 
over the root section and illustrates the progressive change 
to the supersonic type as the Mach number approaches 1.0. 
It can be seen that an increase in Mach number will not only 
increase the distortion of the pressure distribution but will 
increase the extent of the distortion along the span. 



Fiqcf.e 7.— Effect of Macb n umb er on pressure distribution ovor root section of swept-back 
wing, A -RV*, biconvex section. 


An interesting point to be noted is that not all sections of 
the swept-back wing have zero pressure drag. A positive 
drag appears on the root sections and a negative drag on the 
tip sections. Hence the spanwiso drag distribution is 
qualitatively similar to that at supersonic speeds though, of 
course, the net subsonic pressure drag is zero. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., May 1947. 
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APPENDIX 


SYMBOLS 


V flight velocity 

M Mach number 

x, y, z coordinates 
£ point on x axis 

point on y axis 

disturbance-velocity potential 
, v, w disturbance-velocity components 
p local pressure 

g dynamic pressure 

p air density 

P n , Q n Legendre fimctions 
' D differential operator ( dfdx ) 

t thickness of wing 



c 

chord of wing (measured along x) 

m 

slope of line source (absolute value) 

x' 

x+my 

y' 

y—mx 

? 

x — my 

y' 

y+mx 

R. P. 

Real part 
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